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Application of extreme value statistics 
to corrosion is reviewed. It is empha- 
sized that the concept of corrosion 
probability is important for a quantita- 
tive evaluation of corrosion failure and 
its prediction, especially for localized 
corrosion. Extreme value statistics is 
quite useful for assessing the maximum 
pit depth and/or the minimum lime for 
crack generation. The maximum pit 
depth depending on the surface area 
can be evaluated by using ihc Gumbel 
distribution with the concept of return 
period. A standardized procedure is 
proposed for estimating the maximum 
pil depth. Examples of corrosion failure 



analysis using extreme value statistics, 
which were reported mainly in Japan, 
arc shown. Accumulated experiences 
suggest that an appropriate classifica- 
tion of data based on the corrosion 
mechanism is required before applying 
extreme value analysis. 
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1. Introduction 

Development of extreme value statistics and its 
application to various fields, including corrosion, 
has been described by Gumbel in Ref. [1] and [2], 
Evans is one of the pioneers of modern corrosion 
science, and first established the concept of corro- 
sion probability [3, 4], Eldredge [5] used extreme 
value statistics to obtain the maximum value of pit 
depth on an oil well tube wall as a function of tube 
surface area. Scott [6] found a logarithmic depen- 
dence of the maximum pit depth on surface area, 
and explained that dependence by referring to Trip- 
pet [7]. Aziz [8] and Eledredge [5] discussed almost 
all important points to be considered for the analy- 
sis of corrosion pit data and made use of the return 
period concept. This concept, originally introduced 
in the fields of hydrology or meteorology, is now 
used to obtain a size factor which makes it possible 
to estimate the maximum pit depth in a large sur- 
face area based on the distribution of a small num- 
ber of pit depth data from the small surface area 



samples. In Japan, early review papers on corrosion 
probability and extreme value statistics by Masuko 
[9] and Shibata [10, 11] contributed to the study of 
the extreme value statistics as applied to corrosion 
problems. Ishikawa [12,13,14] and Imagawa 
[15,16,17] applied extreme value statistics to ana- 
lyze engineering data. Kase [18, 19] reviewed 
Lieblein's paper [20], introducing MVLUE (mini- 
mum variance linear unbiased estimator) method 
for estimating the distribution parameters, Lieblein 
had given the coefficient of MVLUE up to N = 6. 
Recently, Tsuge [21] had calculated the coefficients 
up to N = 45, and confirmed that the parameters es- 
timated by the MVLUE method are unbiased and 
efficient and are consistent with values estimated by 
the method of moments or maximum likelihood 
when the sample size exceeds more than 20. The 
committee of Japan Society of Corrosion Engineer- 
ing (JSCE) proposed a standard procedure [22] to 
estimate the maximum pit depth from the small 
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sample size data by using the MVLUE method, and 
a computer program named EVAN [23] was devel- 
oped. Recently Laycock et al. [24] reported that a 
generalized extreme value distribution is more con- 
venient for corrosion depth analysis, because no 
preliminary assumption on the type of distribution 
is needed. An introductory book [25] by Kowaka et 
al. helped to differentiate extreme value analysis 
among corrosion workers in Japan. General back- 
ground on extreme value statistics is provided by 
Ang and Tang [26], and Kinnison [27], 



Table 2. Examples of the extreme value analysis for corrosion 
using the Weibull (exponential) distribution 



Example 



Ref. 



Failure life analysis of stress corrosion [56] 

cracking of stainless steel heat exchanger tubes 

Failure life distribution of stainless steels [57] 

in high temperature and high pressure water 

Effect of CaCla concentration on SCC life time [58] 

distributions of stainless steels 

Evaluation of SCC failure life of stainless steel [59] 

in high temperature water 



2. Application of the Extreme Value 
Analysis to Corrosion 

In the early 1980s, meetings and symposia 
[28,29,30,31] were held in Japan for discussing the 
basic principles of extreme value statistics as well as 
difficulties and problems in their application to cor- 
rosion. In Table 1, several topics for which the 
Gumbel distribution is applied are listed. Table 2 
includes cases analyzed using the Weibull distribu- 
tion, including the exponential distribution. Before 
discussing case histories, the standard procedure 
[22] proposed by the committee is briefly explained; 
details are available elsewhere [22, 32]. 

Table 1. Examples of the extreme value analysis for corrosion 
using the Gumbel distribution 



Example 



Ref. 



Life prediction of super heater tubes of the power 

plant 

Application of the extreme value analysis to 

heating tubes of the boiler 

Estimation of the maximum amount of impurity 

segregation in steel 

Failure life estimation of SCC for Ni base alloys 

Extreme value analysis of the corrosion depth of 

the oil tank plate 



10. 

11. 

12. 
13. 

14. 



[49] 

[37] 

[42] 

[41] 
[35] 
[36] 
[50] 

Life prediction of heat exchanger tubes [51] 

Eddy current examination system for heal exchanger [43] 

lubes with the extreme value analysis 

Extreme value analysis of pitting corrosion of heat [52] 

exchanger tubes 

Methods for the parameter estimation of the pit [16] 

distribution in plants 

Ultrasonic method for heat exchanger tubes with [44] 

the extreme value analysis 

Maintenance system for coated heat exchanger tubes [45] 

Corrosion of steels in sea water [53] 

Analysis of perforation of zinc plating steels by [54] 

extreme value statistics 

Fatigue crack behavior of high strength steel in [55] 

artificial sea water 



2,1 The Gumbel Distribution 

The procedure is proposed mainly for analyzing 
pit depth distribution by using the Gumbel distribu- 
tion and the return period in order to estimate the 
maximum depth of the larger surface area from 
which small area specimens are sampled. The Gum- 
bel distribution is expressed as 



F ty) = exp (-exp (- (x-X)ta), 



(1) 



where F(x) is the cumulative probability of pit 
depth, x, and A and a are the location and scale 
parameters. The reduced variate, y, 



y=(jc-A)/a 
is introduced, and then 

y = -ln(-ln(F(»)) 



(2) 



(3) 



is used for constructing the Gumbel probability pa- 
per. Plotting position for the cumulative probability 
can be calculated simply by 



F(y) = l-i/(l+N), 



(4) 



where i is the ith of the ordered value, Jt, in de- 
scending order and N is the total number of sample. 
Plotting y as a function of x yields a best-fitting 
straight line; its slope provides 1/a and its intercept 
(aty =0) yields A. Instead of this graphical estima- 
tion of the parameters, more reliable estimates of a 
and A can be obtained by using the MVLUE (min- 
imum variance unbiased estimator) method, the 
maximum likelihood and the method of moments. 
Among them the MVLUE method which is dis- 
cussed by Lieblein [20] is found to be more efficient 
and unbiased for small size samples. The MVLUE 
estimator can be calculated by 
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X = Xai (N, n)x ; 
a = 2bi (N, n)xi, 



(5) 



where a> (N, n) and bt (N, n) are weights for each 
sample depending on the sample size, N, and trun- 
cated number, n, which are tabulated in the table 
given by Tsuge [21, 23] up to N=45. The weights, 
A, B, C, of variance, V, 

V = a*(A (N,n)y 2 + B(N,n)y+C(N,n)) (6) 

are also found in the table given by Tsuge [21, 23]. 
For the pit depth distribution, the return period, T, 
is defined as 



T=S!s, 



(7) 



where 5 is the surface object (e.g., a tank plate) to 
be examined and s is the area of the small speci- 
mens which are sampled randomly from the objec- 
tive. The return period, T, is in effect a size factor. 
The mode, A, of the pit depth distribution for the 
small specimen is simply obtained by the MVLUE 
estimators mentioned above, and the mode for the 
T times larger surface, x„*x, is given by 



x fflm =A +aln (T). 



(8) 



The perforation probability, P, of the maximum pit 
through the wall thickness, d, is given by 

P = l-exp(-exp(- (d-(k + a\n(T)/a)). (9) 

Finally, the procedure [22] requires reporting the 
surface area of the object, S, with the small sample 
area, s, providing the return period, T ( = S/s), and 
the number of samples, N, with data number, n , ac- 
tually obtained. In addition, the original thickness 
of the plate, d, and the perforation probability, P, 
if needed are to be stated. The above procedure 
does not request to check a goodness of fit of the 
distribution obtained to the Gumbel distribution, 
but recommends to examine the fitness by the Kol- 
mogorov-Smirnov or chi-square test if needed. 

2.2 The Weibull Distribution 

The third type for the smallest value called the 
Weibull distribution 



F(r) = l-exp(-((/-T)Ai)'") 



(10) 



can be fitted to the failure life distribution of stress 
corrosion cracking [33, 34] as shown in Table 2, 



where y, rj and m are the location, scale and shape 
parameter, respectively. This third type of asymp- 
totic distribution for the smallest value can be trans- 
formed to the first type for the largest value, that is, 
Eq. (1), by changing \-F(t) to F{z) and by intro- 
ducing the following reduced variate 



X=]n(t-y),z=(X-\)fa. 



(11) 



The same MVLUE method used for Eq. (1) can be 
utilized for parameter estimation, because the fol- 
lowing relations exist between the parameters of 
both distributions; 



A =ln(i7), a — \lm. 



(12) 



The above unified procedure for estimating 
parameters of the Gumbel and Weibull distribution 
was coded in the computer program EVAN [23]. 



3. Examples 

Several examples are provided to demonstrate 
the usefulness of extreme value statistics for analyz- 
ing corrosion problems. 

3.1 Maximum Pit Depth of Oil Tank Plate 

Through the 1960s and 1970s a number of oil 
tanks were built in Japan. In the late 1970s there 
occurred frequent oil leakages from tanks due to 
corrosion failure. Oil refinery or petrochemical in- 
dustries were located along the seacoast and oil 
leakage caused serious environmental damages. In 
1976, the fire service law was revised to enforce in- 
spection of the thickness of the base and annular 
plates of oil tanks every time oil was evacuated. On 
these occasions extreme value analysis was applied 
and found to be a powerful tool for estimating the 
maximum pit depth. It is emphasized that data for 
the base plate and the annular plate should be con- 
sidered separately because they are characterized 
by different corrosion damage and mechanisms. 

The law requests that plate thickness has to be 
measured at the corners of every 10 cm square on 
the whole surface of the plate. This inspection pro- 
cedure contributed greatly to reducing corrosion 
leakage, but was time-consuming and costly. The 
extreme value analysis was then studied intensively 
in this field [35] [36]. Pit depth distribution sampled 
from the whole base plate was found to obey the 
Poisson distribution. 
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Araki et al. [36] found that the largest value from 
each small square (5 = 1 m 2 ) being randomly sam- 
pled obeys the Gumbel distribution as shown in 
Fig. 1 [36]. The slope and intersect of the line 
(a = 0.694 and A =1.41) were estimated by the 
MVLUE estimators of Eq. (5). In this case, the sur- 
face area of the base plate, S, was 1535 m 2 and the 
return period or size factor was calculated to be 
T = Sis = 1535. The maximum depth, x miX , was cal- 
culated by Eq. (8): 

*™ = 1.41 + 0.694 x In (1535) = 6.50 mm 

which is shown also in Fig. 1. These data were ob- 
tained for a base plate which was exposed for 12.6 
years. Data for both annular and base plate exposed 
for 7.7 years were plotted in Fig. 2 [36], from which 
the first leak due to the maximum pit is to be ex- 
pected after 17.6 years for the base plate and 23.5 
years for the annular plate, respectively. The effect 
of N, s, on the estimates was examined and it was 
concluded that the MVLUE method is optimal for 
N <20 and the maximum likelihood method is reli- 
able for N > 20. 

3.2 Rupture of Heat Exchanger Tubes of the 
Boiler 

Super heater and economizer tubes of boilers are 
exposed to high temperature gases with salt de- 
posits which cause severe corrosion attack. Corro- 
sion attack is not uniform, but localized at several 
sites, wall thinning at the localized site resulting in 
burst. Regular inspection is needed to predict time 
for replacement of the tube before burst. All tubes 
have to be examined for predicting exact time with 
high confidence, but cost of inspection being high, 
that Fukuda et al. [37] introduced the use of 
extreme value analysis to supplement the inspection 
of a small number of tubes. In Fig, 3, the largest val- 
ues of wall thinning observed for 14 tubes are plot- 
ted on Gumbel probability paper. The distribution 
of wall thinning at every inspection time is seen to 
obey the Gumbel distribution and the maximum 
thickness determined by the return period (40 
tubes) increases with operation time as shown in 
Fig. 4. A criterion for a proper replacement time 
has been proposed, which requires replacement 
when the wall thickness reaches half of the design 
thickness, t, y . Risk of burst could be avoided by 
estimating the depth and noting the proposed crite- 
rion. 
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Fig. 1. The Gumbel pint of the maximum pits on (he bottom 
plate of the oil tank. 
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Fig. 2. Estimated depth of the maximum pit on the whole sur- 
face as a function of operation years, and prediction of failure 
life for the bottom and annular plate. 

3.3 The Pit Depth Distribution of Steel Piles In 
Sea Water 

Since the 1970s, steel pipes and piles have been 
used extensively in Japan for harbor construction, 
because lead time for construction could be 
reduced compared with using concrete. Recently, 
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Fig, 3. The Gumbel plots of the maximum thickness loss of 
boiler tubes used for different operation times. 
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Fig. 4. Thickness loss of boiler lubes as a function of operation 
time and estimation of rupture time. 



corrosion of steel pipes and piles was found to cause 
the collapse of harbor structures. Then corrosion 
damage of steel structures exposed in sea water has 
been inspected and analyzed by using extreme value 
analysis. Itoh et al. [38] reported that three differ- 
ent types of depth distribution were found for steel 
piles and plates depending on exposure time and 
exposure location such as water level and deep sea 
(Fig. 5). The type A distribution which exhibits a 
nearly straight line, was found for uniform corro- 
sion loss, its mean value being below 1.0-1.2 mm 
thickness. The type C distribution obeying the 
Gumbel distribution was observed for heavily local- 
ized specimens. The type B distribution is a mixed 
type of A and C distributions. The estimated depth 
using the return period was consistent with observa- 
tions. 
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Fie. 5. Various types of the distribution observed for steel piles 
and pipes exposed in sea water. 

3.4 Classification of Data Based on Corrosion 
Knowledge 

In any of the cases mentioned earlier, measured 
sets of data is fined by two or three distributions 
and must be separated from each other before the 
analysis in order to obtain the maximum value. 
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Imagawa et al. [15,16,17,39] presented many cases 
which require classification of data. For example, 
data for the heat exchanger tubes had to be classi- 
fied into the inlet and outlet side samples because 
corrosion form and its degree of damage were dif- 
ferent at the two locations owing to exposure to dif- 
ferent temperatures. For the oil tank, Imagawa 
observed that more deep pits were formed on the 
welding line compared with other parts. He ob- 
tained the different estimated value of the pit depth 
for each classified sample. At the present time, the 
classification was done on corrosion knowledge and 
experience, but it is required to establish a proce- 
dure based on a common criterion. 

3.5 Crack Depth Distribution of Stress Corrosion 
Cracking 

Stress corrosion cracking is one of the most dan- 
gerous corrosion failure and shows random occur- 



rence which is a very specific and common feature 
of materials fracture. The Weibull distribution has 
been known to be quite useful to analyze the distri- 
bution of fracture strength of various materials [40] 
and also has been found to be applicable for analyz- 
ing failure life distribution due to stress corrosion 
cracking [33, 34]. 

An interesting application of the Gumbel distri- 
bution for analyzing the crack depth distribution 
has been reported by Tsuge [41]. The laboratory ex- 
periment for evaluating the susceptibility of stress 
corrosion cracking of Type 304 stainless steel was 
done by using a bent specimen of u-shape. Bending 
gives stress to the specimen and the environmental 
condition of high pressure water causes many 
cracks, which can be revealed by sectioning the 
specimen after the test as shown in Fig. 6. Distribu- 
tion of the crack depth plotted in the Gumbel prob- 
ability paper showed two lines with an inflection 
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Fig. 6. The distribution of intergntnular corrosion attacks and cracks observed for sensitized 
type 304 stainless steel exposed to the BWR simulated water (DO 8 ppm, 250 °C). 
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point as can be seen in Fig. 7. This inflection point 
was found to correspond just to a depth for initiat- 
ing the intergranular crack. Thus the initiation of 
the intergranular crack growth could be separated 
from the initial process of purely chemical inter- 
granular corrosion. 

3.6 Estimation of the Maximum Segregation of 
Impurities in Steel 

Continuous casting of steel is one of the innova- 
tive technologies achieved by the steel industry. 
Segregation and its band which are formed during 
solidification at the center of slab remain after 
rolling and work as initiation sites for fracture phe- 
nomena such as lamellar tear and hydrogen induced 
cracking (HIC). The maximum amount of segrega- 
tion was found to be related to the above fracture 



phenomena, so that extreme value analysis was ap- 
plied for estimating the maximum amount of segre- 
gation in steel plate from small area samples [42], 
the concentration of impurities being measured by 
using EPMA. The maximum amount of segregation 
thus determined can be used to predict the suscep- 
tibility to lamellar tear. It should be emphasized 
that the statistical procedure for the chemical anal- 
ysis is mainly concerned with the mean and stan- 
dard deviation which assesses the reliability of the 
measurement, but not with extreme values. In re- 
cent years, highly sensitive analytical methods have 
been developed, but it is not clear how to correlate 
the data of small area samples to that of the total or 
bulk specimen. The ratio of the analytical area to 
the bulk specimen reaches almost to 10"*, and ex- 
treme value statistics is expected to be useful. 
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Fig. 7. Distributions changing with exposure time, the initial distribution corre- 
sponding to intergranular corrosion and the second 10 intergranular cracking. 
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3.7 Non-Destructive Methods With the Extreme 
Value Analysis 

Various types of nondestructive methods are 
used for inspecting and examining corrosion dam- 
age in order to prevent failure. High sensitivity and 
resolution in time or in space are required for the 
measurement. In addition, a computer-aided opera- 
tion becomes popular, because huge amounts of 
data must be evaluated. For the heat exchanger, a 
thousand tubes must be checked and the number of 
data easily exceeds 10 s . An eddy current sensor [43] 
and ultrasonic sensor probe [44] to steel tubes, and 
an impedance sensor probe [45] for coating tubes 
have been developed with the data logger and the 
extreme value analysis software. 

4. Discussion 

The size effect on the maximum pit depth is 
found to be estimated with confidence by introduc- 
ing the concept of the return period. Theoretical 
bases of the procedure have been provided by ex- 
treme value theory [1]. Our experience shows that 
the pit depth distribution obeys the normal or expo- 
nential distribution, which belong to the exponen- 
tial distribution family. Thus the maximum values of 
pit depth extracted from the exponential family dis- 
tribution may reasonably be expected to obey the 
Gumbel distribution. Thus the size effect could be 
rationally predicted by using the concept of return 
period. 

Evans [46] pointed out, however, that some cases 
as observed by Wormwell et al. [47] does not obey 
a normal or exponential type distribution, but that 
the tail of the distribution is limited at a certain 
depth. Evans emphasized that such a limit is rea- 
sonable for the case of anodic reaction control situ- 
ation and this limited depth gives a rough indication 
of the greatest pit depth to be expected on a much 
larger area. Evans, however, did not check another 
possibility using the type III distribution which has 
an upper limit. Recently Laycock et al. [24] dis- 
cussed that usefulness of the generalized extreme 
value (GEV) distribution: 



F(jt)-exp (-(l-k(x-u )/«)'* kx < = a+uk, 



(13) 



because the distribution subsumes all three types 
with the sign of a shape parameter, k. When k is 



zero, negative or positive, the distribution changes 
to type I, type II, and type HI, respectively. They 
found that the pit distribution on stainless steels in 
acidified chloride solution fits the GEV distribution 
with k positive, indicating that the type III for the 
largest value could be fitted. The type III distribu- 
tion has a bound or a limit with increasing area, as 
suggested by Evans. 

What sample size, or what size of specimen area 
should be used are questions from non-specialists in 
statistics. For this question, we proposed a proce- 
dure or criterion for choosing s, N and T based on 
the variance given by Eq. (6). The surface area, S, 
of the object is given, and the sampling area, s, is 
selected so as to include at least one pit. Then 
T( = SIs) is obtained. Accumulated data of the 
parameters of a and A suggest [32] that the ratio of 
a/A for localized corrosion is below, or not much 
larger than, 0.3. Kinnison [48] states that the asymp- 
totic theory predicts a constant ratio of 0.313 for all 
extreme value distributions. Then it can be assumed 
that the ratio o/A, is 0.3. If we wish to control vari- 
ance within (A/3) 2 , the following relation can be de- 
duced from Eq. (6) 



(k/a3f=A (N,n)y 1 +B( L N l n)y + C. 



(H) 



Equation (14) can be solved fory or Tas a function 
of N and ff/A, as plotted in Fig. 8. When the ratio 
of a/A can be equated to 0.3 as discussed before, a 
suitable number of samples can be found for a given 
return period, T. From this figure, the required size 
of samples is N = 30 for 7=1200, or N =20 for 
r = 274 and so on. This figure is approximately the 
same as what was observed empirically. 



5. Conclusions 

Extreme value statistics has been found to be a 
powerful tool for estimating the maximum value of 
localized corrosion depending on the surface area. 
Accumulation of data and experience, however, re- 
veals that statistics is less important than corrosion 
experience and knowledge for obtaining a reason- 
able estimation; measured data must be classified 
based on the form of corrosion damage and its de- 
gree before the analysis. Properly classified data is 
found to provide a very reasonable value. Nonde- 
structive methods for measuring wall thickness with 
various types of sensors, combined with extreme 
value analysis, have been developed in recent years. 
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Fig. 8, Calculated curves of the optimum condition tn choose T 
and N at various value of al\. 
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